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Precise measurements of charm production in deep inelastic scattering are pre- 
sented and compared with next-to-leading-order QCD. The data show sensitivity 
to the choice of parametrisation of the gluon density in the proton and could be 
used in future global fits of the parton densities. 



1 Introduction 

Charm quarks are produced copiously in deep inelastic scattering (DIS) at HERA. 
At sufficiently high photon virtualities, Q 2 , the production of charm quarks con- 
stitutes up to 30% of the total cross section |1I2| . Previous measurements of D* 
cross sections [1121314) indicate that the production of charm quarks in DIS in the 
range 1 < Q 2 < 600 GeV 2 is consistent with QCD calculations in which charm is 
produced through the boson-gluon-fusion mechanism. This implies that the charm 
cross section is directly sensitive to the gluon density in the proton. 

In this paper, measurements of the D* cross section are presented with im- 
proved precision and in a kinematic region extending to higher Q 2 than the previous 
ZEUS results p^. The cross sections are compared to a next-to-leading-order (NLO) 
QCD calculation using various parton density functions (PDFs) in the proton. In 
particular, the data are compared to calculations using the recent ZEUS NLO 
fit 0, in which the parton densities in the proton are parametrised by performing 
fits to inclusive DIS measurements from ZEUS and fixed-target experiments. The 
cross-section measurements are used to extract the charm contribution, F^ c , to the 
proton structure function, F^. 



2 Results 

An integrated luminosity of 82 pb _1 was used, more than double that of the previ- 
ous ZEUS result 1 . Events were selected in a DIS regime defined by the kinematic 
region 1.5 < Q 2 < 1000 GeV 2 and the inelasticity, y, restricted to 0.02 < y < 0.7. 
Charm quarks were identified in these events by tagging a D* meson in the decay 
channel D* — > D°ir s — > Kirir s . Transverse momentum, pr(D*), and pseudorapid- 
ity, rj(D*), requirements on the D* meson, governed mainly by the acceptance of 
the tracking detector, were 1.5 < pt{D*) < 15 GeV and \rj(D*)\ < 1.5. After all 
cuts, 5545 ± 129 D* mesons remained for cross-section measurements. 

The D* production rate, r = N/C, in the e~p data set is systematically higher 
than that in the e + p data set. This difference increases with Q 2 ; for example, the 
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ratio of the rates, r e P/r e is equal to 1.12±0.06 for 1.5 < Q 2 < 1000 GeV 2 , while 
for 40 < Q 2 < 1000 GeV 2 it is 1.67 ±0.21 (only statistical errors are given). Such a 
difference in production cross sections is not expected from known physics processes. 
Many checks were performed to try and understand this observed difference; none 
gave an indication of its source. Therefore, the difference in rate is assumed to be a 
statistical fluctuation and the two sets of data were combined for the final results. 
An increased sample of e~p data can resolve the issue. 

The differential D* cross sections as a function of Q 2 , x, pr{D*) and rj(D*) are 
shown in Figure compared with predictions from an NLO QCD calculation 
Predictions are shown for two PDFs: ZEUS NLO and CTEQ5F3 0; and two 
hadronisation schemes: the Peterson function and the Lund string fragmentation [§] 
as in Aroma For all cross sections, the NLO predictions give a reasonable 

description of the data. The NLO calculation does, however, exhibit a somewhat 
different shape, particularly for da/dx, where the NLO is below the data at low x 
and above the data at high x. The predictions using CTEQ5F3 instead of the ZEUS 
NLO fit, or using Aroma for the hadronisation instead of the Peterson function, 
give better agreement with the data. The prediction using the ZEUS NLO fit 
gives a better description than that using CTEQ5F3 (and also better than the 
prediction using GRV98-HO n °t shown) for the cross-section da/drj(D*). A 
better description is also achieved by using Aroma for the hadronisation, although, 
in this case, da/dpx(D*) is not so well described. Further refinement of NLO QCD 
fits and even the use of these data in future fits may achieve a better description. 

The open-charm contribution, i 7 ^, to the proton structure-function F% can be 
defined in terms of the inclusive double-differential cc cross section in x and Q 2 by 

d 2 <j c£ (x,Q 2 ) _ 2ira 2 2 



dxdQ 2 



■{[1 + (1 - yf]Ff (x, Q 2 ) - y 2 F?(x, Q 2 )}. (1) 



In this analysis, the cc cross section is obtained by measuring the D* production 
cross section and employing the hadronisation fraction /(c — > D*) to derive the 
total charm cross section. Since only a limited kinematic region is accessible for the 
measurement of D* mesons, a prescription for extrapolating to the full kinematic 
phase space is needed. The measured F 2 CC in a bin i is given by 



F" (,. Q 2 \ - "'. mcas ( e P Z! D * X ) p' - ,,-.()-) : •)] 

0j,theo(ep — ► L) *X) 

where Ci are the cross sections in bin i in the measured region of pr{D*) and 
r/(D*). The value of F^ theo was calculated from the NLO coefficient functions 0- 
In this calculation, the same parton densities, charm mass (m c = 1.35 GeV), and 
factorisation and renormalisation scales (y^rn 2 + Q 2 ) have been used as for the 
NLO QCD calculation of the differential cross sections. The hadronisation was 
performed using the Peterson fragmentation function. Typical extrapolation factors 
to the full phase space in pr(D*) and rj(D*) are 4.7 at low Q 2 and 1.5 at high Q 2 . 

The data for as a function of x for fixed Q 2 are compared in Figure with 
previous measurements from HI [21 an( i ZEUS TJ and with the ZEUS NLO fit. The 
three sets of data are consistent. The values of F| c are also presented as a function of 
Q 2 at fixed values of x and compared with the ZEUS NLO fit in Figure[21 The data 
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Figure 1. Differential cross sections compared to predictions from NLO QCD. The ratios of the 
cross sections to the central NLO QCD predictions are also shown beneath each plot. 



rise with increasing Q 2 , with the rise becoming steeper at lower x, demonstrating 
the property of scaling violation in charm production. The prediction describes the 
data well for all Q 2 and x. The uncertainty on the theoretical prediction is that 
from the PDF fit propagated from the experimental uncertainties of the fitted data. 
At the lowest Q 2 , the uncertainty in the data is comparable to the PDF uncertainty 
shown. This implies that the double-differential cross sections could be used as an 
additional constraint on the gluon density in the proton. 

3 Conclusions 

The production of D* mesons has been measured in DIS at HERA. Predictions 
from NLO QCD are in reasonable agreement with the measured cross sections, 
which show sensitivity to the choice of PDF and hence the gluon distribution in 
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Figure 2. The measured F£ c as a function of x for fixed Q 2 and as a function of Q 2 for fixed x. 
The current data are compared with previous HI and ZEUS measurements the predictions from 
the ZEUS NLO QCD fit. 



the proton. The ZEUS NLO PDF, which was fit to recent inclusive DIS data, gives 
the best description of the D* data. In particular, this is seen in the cross-section 
da j dr)(D*) . The double-differential cross section in y and Q 2 has been measured 
and used to extract the open-charm contribution to F%. Since, at low Q 2 , the 
uncertainties of the data are comparable to those from the PDF fit, the measured 
differential cross sections should be used in future fits to constrain the gluon density. 

References 

1. ZEUS Coll., J. Breitweg et al, Eur. Phys. J. C12 (2000) 35. 

2. HI Coll., C. Adloff et al, Phys. Lett. B528 (2002) 199. 

3. ZEUS Coll., J. Breitweg et al, Phys. Lett. B407 (1997) 402. 

4. HI Coll., J. Adloff et al, Nucl. Phys. B545 (1999) 21. 

5. ZEUS Coll., S. Chekanov et al, Phys. Rev. D69 (2004) 012004. 

6. ZEUS Coll., S. Chekanov et al, Phys. Rev. D67 (2003) 012007. 

7. B.W. Harris and J. Smith, Phys. Rev. D57 (1998) 2806. 

8. CTEQ Coll., H.L. Lai et al Eur. Phys. J. C12 (2000) 375. 

9. B. Andersson et al Phys. Rep. 97 (1983) 31. 

10. G. Ingelman and J. Rathsman, Comp. Phys. Comm. 101 (1997) 135. 

11. M. Cluck, E. Reya and A. Vogt, Eur. Phys. J. C5 (1998) 461. 



